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Abstract 

We consider a problem of replication of random vectors by ordinary integrals in the 
setting when a underlying random variable is generated by a Wiener process. The goal is to 
find an optimal adapted process such that its cumulative integral at a fixed terminal time 
matches this variable. The optimal process has to be minimal in an integral norm. 
AMS 2000 subject classification: 60H30, 93E20 

Key words and phrases: optimal stochastic control, contingent claim replication, martingale 
representation. 

1 Introduction 

It is well known that random variables generated by a Wiener process can be represented via 
stochastic integrals, as is stated by the classical Clark-Ocone-Haussmann theorem. This result 
leaded to the theory of backward stochastic differential equations and the martingale pricing 
method in Mathematical Finance. 

We consider a problem of replication of random variables by ordinary integrals. The goal 
is to find an optimal adapted process such that its cumulative integral at a fixed terminal time 
matches this variable without error. The optimal process has to be minimal in an integral norm. 
Explicit solution of this problem is found. 
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2 The problem setting and the main result 



Consider a standard probability space (Q,T,P) and standard ti-dimensional Wiener process 
w(t) (with w(0) = 0) which generates the filtration T t = cr{w(r) : < r < t} augmented by all 
the P-null sets in T. For p > 1 and q > 1, we denote by U^ m the class of random processes 
v{t) adapted to Ft with values in R nxm and such that E (||v(')IIl,(o,t)) P < +oo. We denote by 
| • | the Euclidean norm for vectors and Frobenius (i.e., Eucledean) norm for matrices. 

Let / be a /^-measurable random vector, / G L2(Q, Ft, P; R ra ). By the Martingale Repre- 
sentation Theorem, there exists an unique kf G LV;^ such that 

/ = E/+ f k f (t)dw(t). 
Jo 

We assume that there exits 9 G (0, T) such that 

esssupE|/c/(£)| 2 < +oo. 
te[e,T] 

Let g : [0, T) — > R be a given measurable function such there exists c > and a G (0.5, 1) 
such that 

< g(t) < c(T - t) a , git)' 1 <c(l + (T-t)- a ), t€[0,T). (2.1) 

An example of such a function is g(t) = 1 for t < T — r, g(t) = (T — t) a for t > T — r, where 
r G (0, T) can be any number. 

Let U be the set of all processes from L^ 1 such that E J Q g{t)\u{t)\ 2 dt < +oo. 

Let Y{t) be measurable matrix process in R raxra , such that T(t) = g(t)G(t), where G(t) > 
is a symmetric positively defined matrix such that the processes G(t) and G(t) _1 are both 
bounded. Clearly, E J Q T u(t) T T(t)u(t)dt < +oo for u G U. 

Let a G R n , and let b G R raxra be a non-degenerate matrix. 

Consider the problem 

i-t 

Minimize E / u(t) T T(t)u(t) dt over u G U (2.2) 

JO 

%(t) = Ax(t) + bu(t), t€(0,T) 
subject to (2.3) 

ar(0) = a, x{T) = f a.s. 

Note that this problem is a modification of a stochastic control problem with terminal contingent 
claim. These problems were studied intensively in the setting that involve backward stochastic 
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differential equations (BSDE); a first problem of this type was introduced by Dokuchaev and 
Zhou (1999). In this setting, a non-zero diffusion coefficient is presented in the evolution equation 
for the plant process as an auxiliary control process. Our setting is different: a non-zero diffusion 
coefficient is not allowed. Our setting is different since the non-zero diffusion term is not allowed. 
Problem (|2.2p - (j2.3p is a linear quadratic control problem. However, it has a potential to be 
extended on control problems of a general type, similarly to the theory of controlled backward 
stochastic differential equation. 

Lemma 2.1 Let k^(t) = R{t)~ l kf(t) , where 

R{s) = j s Q(t)dt, Q{t) = e A{ - T ~ t hY(t)- l b T e AT{ ?- t \ 
Then kp(-) £ L^ d . 

Theorem 2.2 Problem A2.2\) - (2~3\) has a unique solution in U. This solution is defined as 
u(.) = r(t)- 1 b T e AT( - T ' t ^(t), where a path -wise continuous process such that 

fi(t) = JX + / k^(s)dw(s), 
Jo 

where fi = i?(0) _1 (E/ — q) and q = e AT a. 

Remark 2.3 Restrictions (|2.ip on the choice of T(t) = g(t)G(t) mean that the penalty for the 
large size of u(t) vanishes as t — > T. Thus, we don't exclude fast growing u(t) as t — > T such that 
u(t) is not square integrable. This is why we select the class U of admissible controls to be wider 
than -^22 Xl - I n Dokuchaev (2010), a related result was obtained for a simpler case when it was 
required to ensure that x(T) = E{/|J-e} for some 9 < T. In this setting, the exact match could 
be achieved only for J^-measurable /; the optimal solution was found to be a square integrable 
process. 

3 Proofs 

Proof of Lemma \2.1\ By the assumptions, we have that Q(t) = g{t)~ l Q(t), where 

Q(t)=e A{ - T -%G(ty l b T e AT ^ 
is a bounded matrix, <7(i) -1 > c _1 (T — t)~ a for t 6 [0, T). For a matrix M = M T > 0, set 

p(M) = inf x T Mx. 

x£R": |x|=l 
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We have that £ = inf sg r 0! r] p(Q(s)) > and 



p(R(t)) > f p{Q{s))ds > c" 1 [ T (T - s)- a p(Q(s))ds > -(^—^ 
Jt Jt c 1 - a 



=T 



s=t 



Hence 



It follows from (|3.ip that 



C (T — i) 1 -" 

PW)) < ~ i • (3-1) 

c 1 — a 



(T-t) 1 -"' 

for some constant C > that is defined by c and n. Hence 

£ l^r 1 ! 2 * < ^ 2 (l - «) 2 [ pr^p ^ = C 2 (l - a) 2 (T 2 ~ ^ < +oo. 
It follows that 

E / \k^(t)\ 2 dt<B [ iRit)- 1 ] 2 \k f (t)\ 2 dt 
Jo Jo 

= E iRity^ lkfit^dt + E Ii2(t) -1 | |fc/(t)| 2 cft 
Jo J0 

< sup E / |fc/(i)| 2 di + / |-R(t) -1 | E\k f (t)\ 2 dt 

te[a,e] Jo Je 

r-6 r-T 

= sup |i2(t) _1 | 2 E / (^(t)! 2 ^ + esssupE|A; / (t)| 2 / |22(t) _1 | 2 dt 

tero.01 Jo te\e,T] Je 



and 



E / |^(t)| 2 dt < +oo. (3.2) 
J o 



This completes the proof of Lemma 12.11 □ 

Proof of TheoremlKE By the definition of U, it follows that x(T) G L 2 (0, 7" T , P; R n ) for 
any u € U . Let the function L(u, p) : U x L2(^, Jt, P; -> R be defined as 

L(u, p) = -E / n(t) T r(t)w(t) dt + E/i(/ - x(T)). 
2 Jo 

For a given p, consider the following problem: 

Minimize L(u, p) over it£[J. (3.3) 

We solve problem (|3.3p using the so-called Stochastic Maximum Principle that gives a necessary 
condition of optimality for stochastic control problems (see, e.g., Arkin and Saksonov (1979), 
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Bensoussan (1983), Dokuchaev and Zhou (1999), Haussmann (1986), Kushner (1972)). The only 
u = satisfying these necessary conditions of optimality is defined by 

m = e AT ( T -*)/,(t), fiit) = E{ M |Ji}. (3.4) 

Clearly, the function L(u, /i) is strictly concave in u, and this minimization problem has an 
unique solution. Therefore, this u is the solution of (|3.3|) . 
Further, we consider the following problem: 

Maximize L(u M ,//) over \x G ^(O, Tt-, P; R n ). (3-5) 

Clearly, the corresponding x(T) is 

X (T) = [ T e^-'hu^dt + e AT a, 
Jo 



and 

L(G M , /i) = -E f u M (i) T r(*)v M (t) eft - Bfi T [ e A ^ T -^bu^t)dt - E/x T e AT a + E/i T /. 
2 Jo j 



We have that 

l-T 



Je / u^r&u^dt-EfJ [ e A{T -^bu^t)dt 
* Jo Jo 

= -E f (T(t)- 1 b T i;{t)) T T(t)T{t)- 1 b T i;{t)dt-B^ T [ e A( - T - t ^bT(t)- 1 b T ^j{t)dt 
2 Jo Jo 

= Je/ v(0 T ^r(t)" 1 r(t)r(t)~ 1 6 T ^(t)dt-E^ T / e A(T - t) 6r(t)- 1 6 T ^(t)dt 

2 Jo Jo 

= -e/ ^(t) T bT(ty 1 b T ij(t)dt--Efi T [ e^-HT^y^^dt 
2 Jo Jo 

= -E / T [e AT ( T - t )/i(t)] T &r(t)- 1 6 T e AT ( T - t )^(t) ( it-E / x T [ T e^-'hTity^e^^^dt 
2 Jo Jo 

= ^E l T \(t) T e A ( T -%T(tr 1 b T e Ar ^ T ' t ^(t)dt--Ei, T C \ A ^-%Y(tyW ' e^^^dt 
2 Jo Jo 

= ^E/ /i(t) T g(t)/i(t) dt - E^ T / Q{t)nit)dt = --E n{t) T Q(t)n(t)dt. 
2 Jo Jo 2 Jo 



We used for the last equality that 

rT i-T 



E/x T ! Q(t)fj,(t)dt = E / ii(t) T Q(t)^{t)dt. 
Jo Jo 



It follows that 

L{u p , n) = E/x T (/ - q) - ~E / /i(t) T Q(t)^(t) dt. 

* Jo 

By the Martingale Representation Theorem, there exists k^ G LV£ d such that 

V = fi+ [ k^dwit), (3.6) 



where ji = E/i. We have that 



L(u IM , f i) = fl T {f-q)-lfl T R(0)fl-lE[ dt ( k^(T) T Q(t)k^T)dT + E I k^t) T k f (t)dt 

1 1 Jo Jo Jo 

= fi 1 (f - q) - \^ R(0)fi - J k^TyR(T)k^T)dT + vj^k^t) T k f {t)dt, 



Clearly, the solution of problem (|3.5p is uniquely defined by (|3.6|) with 

fi = R(0y 1 (f-q), k fl (t) = R(ty 1 k f (t). (3.7) 

By Lemma [2TTT it follows from (|3.2I) that sup t6 j ,T] E|/x(t) | 2 < +oo. 
Let show that £ U for any [i. We have that 

( \u»(t)\dt) <CiE( f |r(t) _1 ||/i(t)|dt ) < C 2 sup E|/i(t)| 2 / g{t)~ l dt < +oo. 



E 

In addition 



E 



/" g(t)\ UfM (t)\ 2 dt<C 3 E [ 5 (t)|r(t)-V(0| 2 ^<C4E / 5 (t)-V(i)| 2 * 
o io Jo 

f T 

< C 4 sup E|/i(t)| 2 / g{ty l dt < +oo. 
te[o,T] Jo 

Here Cj > are constants defined by A, b, n, and T. Hence Uu G U. 

We found that sup M inf u L(u, /x) is achieved for (u M ,/i) defined by (|3,6p . (|3.T[) . (|3.4p . We 
have that L(u, /x) is strictly concave in n € C7 and affine in /x G £2(^1 P, R n ). In addition, 
L(u,fi) is continuous in u G L^ 1 given /x G £2^, J 7 , P, R n ), and L(u, fi) is continuous in 
/i G L2(il, J 7 , P, R n ) given By Proposition 2.3 from Ekland and Temam (1999), Chapter 

VI, p. 175, it follows that 

inf sup L(ii,/i) = sup inf L(u, /x). (3-8) 
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Therefore, (u^,//) defined by (|3.6p . (|3.7|) . (|3.4|) is the unique saddle point for 
Let C//be the set of all u(-) £ £/ such that (I2.3P holds. It is easy to see that 

(■T 



1 /" 

inf -E / u(t) T T(t)u(t)dt = inf sup L(u, //), 



and any solution (u, [i) of (|3.8p is such that u 6 It follows that £ Uf and it is the optimal 
solution for problem (|2,2p - (|2,3p . Then the proof of Theorem 12.21 follows. 

4 Example of calculation of u 

Consider a model where / = F(y(T), where y(t) satisfies Ito equation 

dy(t) = a(y(t),t)dt + P{y{t),t)dw{t). 

Here a(x,t) : R n x R — > R n , (3(x,t) : R™ x R — > R nxn are measurable bounded functions such 
that the derivative df5(x,t)/dx is bounded, b(x,t) = ^f3(x,t) (3(x,t) T > 51 > for all x,t, where 

5 > 0, / is the unit matrix. 

Theorem 12.21 can be applied as the following. 

Assume first that a(x, t) = then Theorem 12 . 2 1 ensures that x(T) = f and that E T(t)u(t) 2 dt 
is minimal for 

u(t)=T(t)~ 1 b T e A ( T - t )[ J R(0)- 1 (F(y(0),0)-x(0)) + ^ W" 1 ^ s)dy(s)\ . (4.1) 
Here -ff is the solution of the Cauchy problem for parabolic equation 

^(x,t)+ bij (x, t) f f x (x,t)+J2 a * ( g » f )^ ( g » *) = °» * <T ' 

H(x,T) = F(x). (4.2) 

Here bij,a>i,Xi are the components of the matrix b and vectors /, x. It can be noted that 
H(x,t)=B{F(y(T))\y(t)=x}. 

Assume now that o(-) / 0. We still have that u(t)dt = f for defined by P~T ]) -(|4.2 p : in 
this case, H(x,t) = Eq {F(y(T))\y(t) = x}, where Eq is the expectation is under a probability 
measure Q such that the process y(t) is a martingale under Q. By Girsanov Theorem, this 
measure exists, it is equivalent to the original measure P and unique under our assumptions 
on a and j3. In this case, the value E T(t)u(t) 2 dt is not minimal over u anymore. Instead, 
Eg T(t)u(t) 2 dt is minimal. This still means that deviations of u are minimal but in a different 
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metric. It can be also noted that the definition of the class U for the original measure has to be 
adjusted for the new measure Q, with the expectations E replaced by Eq. 

This model could have practical applications in a number of settings, where it is required to 
ensure that a controlled differentiable process x{t) matches a random vector / generated by a 
uncontrolled stochastic process y(t). For instance, x(t) may represent a controlled path of an 
anti-aircraft missile, and the process y(t) may represent the driving force of an aircraft such that 
/ = F(y(T) represents the aircraft coordinates at time T. 
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